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ABSTRACT 



Magnetic fields are generally expected to increase the characteristic mass of 
stars formed in stellar clusters, because they tend to increase the effective Jeans 
mass. We test this expectation using adaptive mesh refinement (AMR) magneto- 
hydrodynamic simulations of cluster formation in turbulent magnetized clumps 
of molecular clouds, treating stars as accreting sink particles. We find that, con- 
trary to the common expectation, a magnetic field of strength in the observed 
range decreases, rather than increases, the characteristic stellar mass. It (1) 
reduces the number of intermediate-mass stars that are formed through direct 
turbulent compression, because sub-regions of the clump with masses comparable 
to those of stars are typically magnetically subcritical and cannot be compressed 
directly into collapse, and (2) increases the number of low-mass stars that are 
produced from the fragmentation of dense filaments. The filaments result from 
mass accumulation along the field lines. In order to become magnetically su- 
percritical and fragment, the filament must accumulate a large enough column 
density (proportional to the field strength), which yields a high volume density 
(and thus a small thermal Jeans mass) that is conducive to forming low-mass 
stars. We find, in addition, that the characteristic stellar mass is reduced further 
by outflow feedback. The conclusion is that both magnetic fields and outflow 
feedback are important in shaping the stellar initial mass function (IMF). 

Subject headings: ISM: magnetic fields — Magnetohydrodynamics (MHD) - 
ISM: clouds 
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1. Introduction 

How the initial mass function (IMF) of stars originates is one of the most basic questions 
that a complete theory of star formation must answer. It is also one of the most difficult 
(see reviews by Bonnell et al. 2007 and McKee & Ostriker 2007 and references therein). 
Many ideas have been advanced to explain the IMF, including Zinnecker (1982), Adams 
& Fatuzzo (1996), Elmegreen (1997), Padoan & Nordlund (2002), Larson (2003), Shu et 
al. (2004), Hennebelle & Chabrier (2008), and Kunz & Mouschovias (2009), among others. 
Particularly intriguing is the proposal by Padoan & Nordlund (2002) that the IMF is essen- 
tially determined by the supersonic turbulence, which controls the density distribution in 
the highly inhomogeneous clouds. They used the Jeans criterion to determine analytically 
the mass distribution of the gravitationally unstable regions, which was taken to represent 
the stellar IMF. Both ingredients of the theory, the turbulent structuring of cloud density 
and the criterion for gravitational collapse, are strongly affected, however, by a dynamically 
important magnetic field. Some magnetic effects on the IMF, such as the magnetic cushion 
of turbulent compression (Padoan & Nordlund 2002) and the additional cloud support by 
magnetic pressure (Hennebelle & Chabrier 2008), can be incorporated into the analytic the- 
ory approximately. Other equally, if not more, important magnetic effects are not captured 
by the theory, and will be the focus of our investigation. 

Our study will concentrate on the parsec-scale, cluster-forming dense clumps of molec- 
ular clouds, where the majority of stars, especially massive stars, are thought to originate 
(Lada & Lada 2003), and where magnetic fields have been observed in some cases. For ex- 
ample, in the nearest region of active massive star formation, OMC1, a well-ordered field is 
inferred from polarized dust emission (e.g., Vaillancourt et al. 2008). Crutcher et al. (1999) 
obtained, from CN Zeeman measurements, a line-of-sight field strength Bi os = 360 /iG for 
this region, which corresponds to a dimensionless mass-to-flux ratio A = 2irG 1 / 2 M/& ~ 4.5. 
The ratio is close to the median value A ~ 6 obtained by Falgarone et al. (2008) for a 
sample of dense clumps of massive star formation. Correcting for the projection effects may 
reduce the ratio by a factor of 2-3 (Shu et al. 1999), yielding a median value for the intrinsic 
mass-to-flux ratio of 2-3. 

A magnetic field corresponding to A ~a few would not be strong enough to support the 
clump against gravitational collapse by itself. It can, however, change the mass distribution 
of the stars formed. This is because individual stars form out of sub-regions of the clump, 
which are generally more strongly magnetized (relative to their masses) than the clump as 
a whole (i.e., with smaller values of A, Tilley & Pudritz 2007; Dib et al. 2007), because the 
mass of a region (which is proportional to volume) decreases faster with its size than the 
magnetic flux (which is proportional to area). In particular, sub- regions smaller than the 
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clump by a factor comparable to A are expected to be magnetically subcritical (with A < 1) 
in general. It would be hard to directly compress such regions into collapse by turbulence in 
the ideal MHD limit (see, however, Nakamura & Li 2005 and Basu et al. 2009 for studies that 
include ambipolar diffusion). Another effect is that the magnetic forces act on the turbulent 
flows anisotropically, channeling matter along the field lines into dense, flattened structures 
that can subsequently fragment thermally. 

In this Letter, we seek to quantify the above effects using the cluster formation simula- 
tions of Wang et al. (2010), which were carried out with an AMR MHD code that includes 
sink particles and outflow feedback. We find, in § [2], that the characteristic mass of the IMF 
is lowered by both the magnetic field and outflow feedback. The lowering of the characteris- 
tic stellar mass by magnetic fields is somewhat counter-intuitive. It is interpreted physically 
using a so-called "magnetically critical Jeans mass" (eq. [0]) in § El 

2. Result 

In Wang et al. (2010), we simulated cluster formation in moderately condensed, isother- 
mal (T = 20 K), clump of 2 pc in size and 1641 M in mass. We added to the clump, one 
by one, an initial turbulent velocity field of rms Mach number 9 (Model HD), an initially 
uniform magnetic field of 100 /iG (corresponding to a dimensionless mass-to- flux ratio of 1.4 
for the clump as a whole and a larger value of 3.3 for the central, high-density part; Model 
MHD), and mechanical feedback from collimated outflows driven by accreting stellar objects 
(Model WIND). For simulation details we refer the reader to Wang et al. (2010), where we 
also presented results on the global star formation rate and especially the formation of the 
most massive object in each simulation. Here, we will concentrate on the distribution of the 
stellar masses, with an emphasis on low and intermediate mass objects. 

Figure 1 shows the stellar mass distributions for all three models. It is a log-linear plot 
of the number of stars in logarithmic mass bins of A[log(M)] = 0.2 (where M is the stellar 
mass in M ). The distributions are all evaluated at the time when ~ 16% of the total mass 
has been converted into stars. It corresponds to 1.16 Myr for Model WIND (the last output), 
and 0.576 and 0.723 Myr for Model HD and MHD, respectively; for comparison, the free-fall 
time is 0.210 Myr at the clump center initially and 0.475 Myr in a sphere of 1 pc in radius. 

Before proceeding to analyze the numerically obtained mass distributions in Fig. 1, we 
should caution the reader that the mass of a star can be affected by many factors that are 
only crudely modeled in the current generation of large-scale cluster formation simulations, 
including our own. These include (1) outflow feedback, which can potentially remove a 
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Fig. 1. — The number of stars in logarithmic mass bins of A[log(M)] = 0.2 for Model HD 
(dashed), MHD (solid) and WIND (dotted), when 16% of the clump mass has been converted 
into stars. 
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large fraction, perhaps the majority of the mass of a star-forming core (Shu et al. 1987; 
Matzner & McKee 2000; Myers 2008), (2) magnetic fields, which affect not only how dense 
cores form but also how they collapse, and (3) radiative feedback which, once a protostellar 
system (protostar plus disk) appears, heats up the circumstellar material and suppresses 
fragmentation close to the central object (e.g., Krumholz et al. 2007; Offner et al. 2009; 
Price & Bate 2009; Smith et al. 2009; Urban et al. 2009). It is currently not feasible to 
treat simultaneously all these factors in detail. Our strategy is to focus on magnetic fields 
and outflow feedback, which are less well studied, and treat the effects of radiative feedback 
as simply as possible. 

Specifically, we crudely capture the suppression of fragmentation near a star due to 
radiative feedback using a sink-particle merging algorithm, following Krumholz et al. (2004, 
see also discussion in Federrath et al. 2010). It eliminates very low mass (< M mcrg = 
0.01 M ) sink particles within a distance Z mcr g = 1000 AU of an existing star (with mass 
> M merg ). Experimentation shows that increasing M merg to 0.1 M does not change the star 
formation rate or the stellar mass distribution significantly. The adopted merging distance 
^merg is comparable to the radius of the so-called "sphere of thermal influence" (where the 
temperature of the radiatively heated protostellar envelope drops to the ambient value) of 
an object of solar mass and luminosity (see equation [14] and Figure 2 of Adams & Shu 
1985). We find that increasing Z merg from 1000 AU to 2000 AU has little effect on the star 
formation rate and stellar mass distribution. Decreasing it to 600 AU (or three times the 
size of the finest cell) does not change the total star formation rate, but can change the 
number and masses of stars significantly, by up to 50%. We believe, however, that the 
smaller merging distance is less realistic, since the gas on the smaller scale is expected to be 
more strongly heated. The larger merging distance prevents a collapsing core from breaking 
up into many small pieces, in agreement with the fact that the best observed cores in nearby 
star-forming regions typically harbor one, at most a few, stellar systems (e.g., Mundy et al. 
2000). Nevertheless, in view of our crude representation of radiative feedback, we elect to 
focus on the difference between the stellar mass distributions of the three simulations, which 
use the same sink particle treatment, rather than the distributions themselves. 

The magnetic field is solely responsible for the difference between Model HD and MHD. 
From Fig. 1, we see that its main effects are to reduce the number of those intermediate-mass 
stars with masses around ~ 4 M , and to increase the number of lower mass stars with masses 
around ~ 1 M . To be more quantitative, we fit the high mass end of the distribution with a 
power-law, dN/dlog(M) oc M r , as shown in Fig. [2J In all three cases the mass distribution 
deviates sharply from the power-law fit below a characteristic mass, M c ^ (sometimes dubbed 
the "knee" of the IMF, e.g., Bonnell et al. 2006). The values of M c h and power-law index T 
are listed in Table 1. Note that the characteristic mass of ~ 5 M in Model HD is not far 
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from the initial Jeans mass near the clump center, ~ 9.5 M , consistent with Bonnell et al. 
(2006). The power-indexes are not well determined because of the limited number of stars 
in the power-law regime. Nevertheless, they are all comparable to the Salpeter value within 
uncertainties. The characteristic masses are, however, significantly different. In particular, 
the magnetic field in Model MHD has apparently lowered the characteristic mass by a factor 
of ~ 4 compared to that in Model HD. 

To help understand why the field hampers the formation of intermediate-mass stars, we 
show in Figure [3] a snapshot of the so-called "HD excess stars," defined somewhat arbitrarily 
as those stars in HD model in Figs. 1 and 2 with masses in the three mass bins between 
10 1 ' 2 and 10 1 ' 8 M . The majority of these stars are produced at relatively isolated regions, 
by localized converging flows in the initial turbulence. Such regions are largely absent in 
the MHD model, where a moderately strong global magnetic field is present. The reduction 
in the number of intermediate-mass stars makes the stellar mass distribution significantly 
narrower in the MHD model than in the HD model (see Table 1 and Figs. 1 and 2). 

The magnetic field in the MHD model interacts strongly with the turbulence. It changes 
the "texture" of the turbulent clump, by not only resisting the cross-field turbulent com- 
pression, but also promoting mass accumulation along the field lines into flattened sheets or 
filaments. Fig. H] shows an example of the (flattened) filaments resulted from the intrinsically 
anisotropic magnetic forces. It is in such magnetically-induced filaments that the majority 
of the stars in the MHD model form. As we discuss in § [31 the high density in the fila- 
ments favors the formation of low-mass stars, which are more numerous compared to the 
non-magnetic HD model. 

Outflow feedback is responsible for the difference between Model WIND and MHD. 
From Figs. 1 and 2, it is clear that the outflows have (1) shifted the characteristic mass 
M c h to a lower value, from ~ 1.3 M to ~ 0.5 M , and (2) nearly doubled the number of 
stars, from 119 to 203 ( see Table 1). The lowering of the stellar mass by outflow feedback is 
perhaps not too surprising, because, individually, the outflow from a given star can remove 
part of the dense envelope that feeds the growing star (Shu et al. 1987; Matzner & McKee 
2000; Myers 2008) and, collectively, the feedback from all stars can prevent the clump from 
rapid global collapse (Nakamura & Li 2007; Wang et al. 2010), lowering the total rate of 
stellar mass accretion. There is, however, a large spread in both the rate and the duration 
of the mass accretion in stars of similar masses. The spread makes it difficult to develop 
a complete picture of how the outflows lower the characteristic stellar mass precisely. The 
number of stars produced per unit time (i.e., the stellar production rate) turns out to be 
insensitive to the outflow feedback. It takes longer, however, to convert the same amount of 
gas into stars in the presence of the feedback, which is the main reason for the larger number 
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Fig. 2. — Power- law fits to the high mass end of the stellar mass distributions of the three 
models: HD (dashed lines), MHD (solid) and WIND (dotted). The top (bottom) curve is 
raised (lowered) by a factor of 20 for clarity. 
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Fig. 4. — 3D view of a flattened filament in the MHD model. Plotted are isodensity surfaces 
and magnetic field lines (white) at a relatively early time 0.262 Myr. The high density of 
the filament favors the formation of low-mass stars. 
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of stars in Model Wind than in Model MHD. 



3. Discussion and Conclusion 

We have shown that a magnetic field of strength in the observed range changes the IMF 
significantly, by impeding direct turbulent compression of relatively low density material into 
relatively massive stars and, more importantly, promoting the formation of dense, flattened 
structures that subsequently fragment preferentially into low-mass stars. 

The basic reason for the magnetic field to impede direct turbulent fragmentation is that 
a sub-region of the clump that is magnetically subcritical cannot be compressed into prompt 
collapse by turbulence. For a region of size / in a clump of average field strength B and 
mass density po, the characteristic dimensionless mass-to- flux ratio is 

2ttG 1 / 2 M 1 27iG l ^p l 3 2txG 1 ' 2 Po I . . 

Xl = — l 5^ = ^^' W 

which indicates that regions smaller than the characteristic size 

are typically magnetically subcritical (with A < 1). The characteristic size is related to the 
clump size L by Ib ~ L/X cp (where "cp" stands for "clump"), since the mass-to-flux ratio 
of the clump as a whole is A cp ~ 2tcG 1 ^ 2 poL/B , according to equation (CQ). In other words, 
a typical region that is smaller than the clump by a factor of ~ A cp cannot be compressed 
directly into collapse by turbulence. The corresponding characteristic mass is 

where we have scaled the clump mass M cp by 10 3 M Q , and the clump mass-to- flux ratio 
by 2.5, comparable to the (projection-corrected) median value inferred from CN Zeeman 
observations (see discussion in §[[]). 

The characteristic mass Mb is essentially the magnetic Jeans mass at the average clump 
density (see equation [41] of McKee 1999, which gives a numerical value three times larger). 
It is to be compared with the average thermal Jeans mass 
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(e.g., Spitzer 1978). The average magnetic Jeans mass Mb is therefore typically larger than 
the average thermal Jeans mass Mj, indicating that the formation of those stars with masses 
comparable to Mj or smaller by direct turbulent compression can be strongly hindered by 
the magnetic field. The most massive stars are expected to be least affected by the magnetic 
field. However, their formation may be governed more by the global clump collapse (retarded 
by outflow feedback and magnetic fields) than by direct turbulent compression (Wang et al. 
2010; see also Smith et al. 2009). Indeed, massive (> 1OM ) stars in Model MHD contain 
more mass than those in Model HD, indicating that their formation is not hampered by 
magnetic fields, unlike intermediate mass stars. 

The magnetic field can also promote fragmentation, because of the magnetic forces are 
intrinsically anisotropic. It is easier for the turbulent flows to move material along the field 
lines than across them. The net result is the formation of flattened, sheet- or filament-like 
condensations, which are conducive to fragmentation. Larson (1985) estimated the Jeans 
mass for a self-gravitating sheet (thermally supported in the vertical direction) to be 

4 -67 a 4 

M JiSh « (5) 

where a is the isothermal sound speed, and E the column density. In order to form stars, 
the sheet must first become magnetically supercritical, with a column density greater than 
the critical value S cr = _B /(27rG 1 / 2 ). Substituting S cr into equation we have 

9.34 vr a 4 / L \ 2 /10 3 M \ ( T \ 2 /A cp 



M,„ « = 0.78 ^— J ^ (Me) (•) 

which is the mass expected for stars formed through thermal fragmentation of magnetically 
critical filaments. A similar, but somewhat smaller, characteristic mass M = 7r 2 a 4 / (G 3 ^ 2 B ) 
was defined in Shu et al. (2004) in their study of the collapse of magnetized singular isother- 
mal spheres. The important conceptual point is that the magnetically critical Jeans mass 
Mj,cr (or Mo) decreases with an increasing magnetic field strength, because a stronger mag- 
netic field can support a higher column density, which in turn leads to a higher volume 
density (due to gravitational compression along the field lines) and thus a lower thermal 
Jeans mass. It should be relatively insensitive to turbulence, which is expected to be weak 
at the relevant (high) densities. In any case, the magnetically critical Jeans mass is typically 
smaller than the average thermal Jeans mass, by a factor 



M*„ 0ia (JLY(^)"(—Y'(±\ (7) 
Mj ' \lpcj V A/ cp ) \wk) \2.5j k ' 

which is relatively insensitive to clump parameters. It is therefore not too surprising that 
a moderately strong field can lower the characteristic mass of the stars formed in a cluster 
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significantly, as we found numerically. Magnetism is another way of setting a stellar mass 
scale that is distinct from the one advocated by Larson (2005), based on detailed thermal 
physics of molecular cloud material (e.g., Jappsen et al. 2005 and Bonnell et al. 2006). 

To summarize, we have demonstrated through both AMR MHD simulations (§ |2]) and 
analytic considerations (eqs. [l|-[Tj) that moderately strong magnetic fields of the observed 
strengths (corresponding to a dimensionless mass-to-flux ratio of a few) are important in 
shaping the mass distribution of stars produced in turbulent, cluster-forming, dense clumps. 
The basic reason is that the field interacts strongly with the turbulent flows, which affects 
how the flows create regions unstable to gravitational collapse and star formation. The 
magnetic field prevents sub-regions of the clump that are magnetically subcritical (with 
masses less than the average magnetic Jeans mass Mb) from being compressed directly into 
collapse by the turbulence. The magnetic resistance to turbulent compression has apparently 
reduced the number of intermediate mass stars formed in our simulations. More importantly, 
the intrinsically anisotropic magnetic forces channel part of the clump material along the 
field lines into dense structures, where the high density promotes the formation of low-mass 
stars with masses comparable to the magnetically critical Jeans mass Mj cr , which lowers the 
characteristic stellar mass. The mass is further reduced by outflow feedback, which affects 
not only the individual stars that drive the outflows, but also the dynamics of the clump as 
a whole. We conclude that magnetic fields and outflow feedback are important factors that 
should be accounted for in a complete theory of the stellar IMF. 

This work was supported in part by NASA grants (NNG06GJ33G and NNX10AH30G) 
and a Grant-in- Aid for Scientific Research of Japan (20540228). 



Table 1. Parameters of Stellar Mass Distribution 



Model 


IMF exponent T 


Characteristic mass M c h 




Star number 


HD 


-1.2 


5.0 M© 


0.8 


96 


MHD 


-1.4 


1.3 M & 


0.8 


119 


WIND 


-1.0 


0.5 M 


1.8 


203 



Note. — a) Full width of log(M) at half maximum of the stellar mass distribution. 
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